B样条方法在表示与设计自由型曲线曲面形状时显示了强大的威力,然而 | |
在表示与设计初等曲线曲面时时却遇到了麻烦。因为B样条曲线包括其特例的 | |
Bezier曲线都不能精确表示出抛物线外的二次曲线,B样条曲面包括其特例的 | |
Bezier曲面都不能精确表示出抛物面外的二次曲面,而只能给出近似表示。 | |
提出NURBS方法,即非均匀有理B样条方法主要是为了找到与描述自由型曲线 | |
曲面的B样条方法既相统一、又能精确表示二次曲线弧与二次曲面的数学方法。 | |
NURBS方法的主要优点: | |
(1)既为标准解析形状(即前面提到的初等曲线曲面),又为自由型曲线 | |
曲面的精确表示与设计提供了一个公共的数学形式。 | |
(2)修改控制顶点和权因子,为各种形状设计提供了充分的灵活性。 | |
(3)具有明显的几何解释和强有力的几何配套技术(包括节点插入、细 | |
分、升阶等)。 | |
(4)对几何变换和投影变换具有不变性。 | |
(5)非有理B样条、有理与非有理Bezier方法是其特例。 | |
不过,目前应用NURBS中还有一些难以解决的问题: | |
(1)比传统的曲线曲面定义方法需要更多的存储空间,如空间圆需7个参 | |
数(圆心、半径、法矢),而NURBS定义空间圆需38个参数。 | |
(2)权因子选择不当会引起畸变。 | |
(3)对搭接、重叠形状的处理很麻烦。 | |
(4)反求曲线曲面上点的参数值的算法,存在数值不稳定问题。 | |
3.4.1 NURBS曲线的定义 | |
NURBS曲线是由分段有理B样条多项式基函数定义的: | |
|
|
其中,Ri,k(t)(i=0,1,…,n)称为k阶有理基函数,Ni,k(t)是k 阶B样条 | |
基函数,Pi(i=0,1,…,n)是特征多边形控制顶点位置矢量;w i是与Pi对应 | |
的权因子,首末权因子w 0,w n>0,其余w i3 0,以防止分母为零及保留凸包 | |
性质、曲线不因权因子而退化为一点;节点矢量为T=[t0, t1, … , ti, …, | |
点tn+k],节个数是m=n+k+1(n为控制项的点数,k为B样条基函数的阶数)。 | |
对于非周期NURBS曲线,常取两端节点的重复度为k,即 | |
有:,在大多数实际应用中,a =0, | |
b =1。P(t)在区间上是一个k-1次有理多项式,P(t)在整条曲线上 | |
具有k-2阶连续性,对于三次B样条基函数,具有C2连续性。当n=k-1时,k阶 | |
NURBS曲线变成k-1次有理Bezier曲线,k阶NURBS曲线的节点矢量中两端节点的 | |
成节点重复度取k+1就使得曲线具有同次有理Bezier曲线的端点几何性质。 | |
Ri,k(t)具有k阶B样条基函数类似的性质: | |
(1)局部支承性:Ri,k(t)=0,t? [ti, ti+k]; | |
(2)权性:; | |
(3)可微性:如果分母不为零,在节点区间内是无限次连续可微的,在 | |
节点处 (k-1-r)次连续可导,r是该节点的重复度。 | |
(4)若w i=0,则Ri,k(t)=0; | |
(5)若w i=+¥ ,则Ri,k(t)=1; | |
(6)若w j=+¥ ,且j1 i,则Ri,k(t)=0; | |
(7)若w j=1,j=0,1,…,n, 则是B样条基函数;若w | |
jj=1,=0,1,…,n,且 则,Bi,k(t)是 | |
Bernstein基函数。 | |
Ri,k(t)与Ni,k(t)具有类似的性质,导致NURBS曲线与B样条曲线也具有类 | |
似的几何性质: | |
(1)局部性质。k阶NURBS曲线上参数为的一点 | |
至多与k个控制顶点Pi及权因子有关,与其它顶 | |
点和权因子无关;另一方面,若移动k次NURBS曲线的一个控制顶点Pi或改变所 | |
联系的权因子仅仅影响定义在区间上那部分曲线的形状 | |
(2)变差减小性质。 | |
(3)凸包性。定义在非零节点区间上曲线段 | |
位于定义它的k+1个控制顶点的凸包内。整条NURBS曲线位于所 | |
有定义各曲线段的控制顶点的凸包的并集内。所有权因子的非负性,保证了 | |
凸包性质的成立。 | |
(4)在仿射与透射变换下的不变性。 | |
(5)在曲线定义域内有与有理基函数同样的可微性。 | |
(6)如果某个权因子为零,那么相应控制顶点对曲线没有影响。 | |
(7)若,则当时,。 | |
(8)非有理与有理Bezier曲线和非有理B样条曲线是NURBS曲线的特殊情 | |
况 。 | |
3.4.2 齐次坐标表示 | |
为了便于讨论,我们考虑平面NURBS曲线的情况。图3.1.34所示,如果给 | |
一组控制顶点及对应的权因子 | |
则在齐次坐标系xyw中的控制顶点为。 | |
齐次坐标下的k阶非有理B样条曲线可表示为: | |
|
|
若以坐标原点为投影中心,则得到平面曲线: | |
|
|
|
|
三维空间的NURBS曲线可以类似地定义。即对于给定的一组控制顶点 | |
及对应的权因子,则有相应 | |
带权控制点,定义了一条四维的 | |
k阶非有理B样条曲线,然后,取它在第四坐标的超平面上的中心 | |
投影,即得三维空间里定义的一条k阶NURBS曲线。这不仅包含了明确的 | |
几何意义,也说明,非有理B样条的算法可以推广到NURBS曲线,只不过是在 | |
齐次坐标下进行。 | |
3.4.3 权因子的几何意义 | |
由于NURBS曲线权因子w i只影响参数区间定义在区间 | |
上的那部分曲线的形状,因此,我们只考察整条曲线的这一部分。如果固定曲线的参数t,而使变化,则NURBS曲线方程变成以为参数的直线方程,即 | |
NURBS曲线上t值相同的点都位于同一直线上,如图3.1.35所示。我们把曲线与 | |
有理基函数的记号用用如下包含其权因子为变量的记号替代。因当时 | |
,的,故该直线通过控制顶点, | |
分别是对应曲线上的点,即, | |
,。 | |
令 a =Ri,k(t; w i=1 ),b = Ri,k(u) |
|
N,Bi可表示为: | |
|
|
用a 、b 可得到下述比例关系: | |
|
|
上式是(Pi,Bi,N,B)四点的交比,由此式可知: | |
(1)若w i增大活减小,则b 也增大或减小,所以曲线被拉向或推离开 | |
Pi点; | |
(2)若w j增大或减小,曲线被推离或拉向Pj(j1 i)。 | |
|
|
3.4.4 圆锥曲线的NURBS表示 | |
若取节点向量为,则NURBS曲线退化为二次Bezier曲线, | |
且 | |
可以证明,这是圆锥曲线弧方程,称为形状因子,的值 | |
确定了圆锥曲线的类型。时,上式是抛物线弧,时, | |
上式是双曲线弧,时,上式是椭圆弧。且时,上式 | |
退化为一对直线段和,时,上式退化为连接 | |
两点的直线段,如图3.1.36所示。 | |
|
|
3.4.5 NURBS曲线的修改 | |
NURBS曲线的修改有多种方式,常用的方法有修改权因子、控制点和反插 | |
节点。 | |
1.修改权因子 | |
权因子的作用是:当保持控制顶点和其它权因子不变,减少或增加某权因 | |
子时,曲线被推离或拉向相应顶点。假定已给k阶(k-1)次NURBS曲线上参数 | |
为t的一点S,欲将曲线在该点拉向或推离控制顶点一个距离d,以得到新点 | |
,可由重新确定相应的权因子使之改变为来达到,如图3.1.37所示。 | |
|
|
|
|
其中,表示和两点间的距离,d有正负之分,若在和之 | |
间,即曲线被拉向顶点和,d为正,反之为负。 | |
修改过程是拾取曲线上一点,并确定该点的参数,再拾取 | |
控制多边形的一个顶点,它是k+1个控制顶点中的一个,即 | |
,便可算出两点间的距离d。若在直线段上拾取一个点 | |
,就能确定替代老权因子的新权因子,修改后的曲线将通过点。 | |
2.修改控制顶点 | |
若给定曲线上参数为的一点S,方向矢量V和距离d,计算控制顶点的 | |
新位置,以使曲线上S点沿V移动距离d到新位置。可表示为: | |
|
|
于是: |
|
由此可得新控制顶点: |
|
3.反插节点 | |
给定控制多边形顶点与权因子及节点矢量 | |
,就定义了一条k阶NURBS曲线。现欲在该多边形的 | |
的边上选取一点,使得点成为一个新的控制顶点,这就是所谓反插节点。 | |
点可按有理线性插值给出: | |
|
|
于是: |
|
所以 |
|
这就是使得成为一个新控制顶点而要插入的新节点。 | |
当插入新节点使成为新控制顶点的同时,将有k-2个老控 | |
制顶点被包括在内的新控制顶点所替代,如图3.1.38所示。 | |
|
|
3.4.6 非均匀有理B样条(NURBS)曲面 | |
1.NURBS曲面的定义 | |
由双参数变量分段有理多项式定义的NURBS曲面是: | |
|
|
式中是矩形域上特征网格控制点列,是相应控制点的权因子,规定 | |
四角点处用正权因子,即,其余。 | |
和是p阶和q阶的B样条基函数,是双变量有理基函数: | |
|
|
节点矢量和按de Boor递推 | |
公式决定,通常具有下面的形式: | |
p个 q个 p个 q个 |
|
2.NURBS曲面的性质 | |
有理双变量基函数与非有理B样条基函数相类似的性质: | |
(1)局部支承性质:,当或 | |
; | |
(2)权性:; | |
(3)可微性:在每个子矩形域内所有偏导数存在,在重复度为r的u节点 | |
处沿u向是p-r-1次连续可微,在重复度为r的v节点处沿v向是q-r-1次连续可 | |
微; | |
(4)极值:若p,q>1,恒有一个极大值存在; | |
(5)是双变量B样条基函数的推广。 | |
NURBS曲面与非有理B样条曲面也有相类似的几何性质,权因子的几何意义 | |
及修改、控制顶点的修改等也与NURBS曲线类似,这里不在赘述。 | |
我们已经知道,计算机中表示形体,通常用线框、表面和实体三种模 | |
型。线框模型和表面模型保存的三维形体信息都不完整,只有实体模型才能 | |
够完整地、无歧义地表示三维形体。前面我们已经介绍了曲线曲面常用的的 | |
表示形式及其理论基础,从本小节开始,我们介绍实体造型技术的有关问题, | |
主要包括形体在计算机内的表示、分类求交算法和典型的实体造型系统。 | |